Unexpected features of quantum degeneracies in a pairing model with two integrable 

limits 
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The evolution pattern of level crossings and exceptional points is studied in a non-integrable pair- 
ing model with two different integrable limits. One of the integrable limits has two independent 
parameter-dependent integrals of motion. We demonstrate, and illustrate in our model, that quan- 
tum integrability of a system with more than one parameter-dependent integral of motion is always 
signaled by level crossings of a complex-extended Hamiltonian. We also find that integrability im- 
plies a reduced number of exceptional points. Both properties could uniquely characterize quantum 
integrability in small Hilbert spaces. 
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The search for fingerprints of the chaotic/regular dy- 
namics in the quantum regime is often focused on studies 
of spectral properties of the quantum systems. In this 
context, spectral fluctuations were intensely studied for 
various quantum systems. These studies lead to the BGS 
conjecture [l[ that in the semiclassical limit the spectral 
fluctuations of chaotic systems are described by random 
matrix theory. For quantum integrable systems, Berry 
and Tabor [2|] showed that the spectral fluctuations are 
well described by a Poisson statistic. While chaotic sys- 
tems are characterized by level repulsion between succes- 
sive levels, levels of integrable systems are uncorrelated 
allowing crossings between states of the same symme- 
tries. 

Level crossings and degeneracies are important for the 
understanding of spectral fluctuations Q and the onset 
of quantum chaos [4|- Much effort has been devoted to 
studies of degeneracies associated with avoided crossings 
in quantal spectra, focusing mainly on the topological 
structure of the Hilbert space and the geometric phases 
[H, Among these degeneracies, one finds a diabolic 
point where two Riemann sheets of eigenvalues touch 
each other [H, 0] , and an exceptional point (EP) [1, [tj [1(3] 
where the two sheets are entangled by the square-root 
type of singularity. EP appears in the complex g plane of 
a generic Hamiltonian H(g) = Hq + gHi, where both Ho 
and Hi are hermitian and [H ,Hi] ^= 0. In many-body 
systems, EPs have been studied in schematic models like 
the Lipkin model and the interacting boson model 
along the line connecting the dynamical symmetries U(5) 
to 0(6) [HI- Both models can be considered as particu- 
lar examples of 2-level boson pairing models pertaining to 
the general class of integrable Richardson-Gaudin (RG) 
models 13, llJ] . In spite of the fact that level crossings of 



eigenstates with the same global symmetries are only al- 
lowed for integrable systems, 2-level pairing models have 
no level crossing and all degeneracies take place in the 
complex plane as EPs. 

In this work, we will introduce a prototypical quan- 
tum integrable system, the 3-level RG model, to dis- 
cuss the appearance of level crossings and EPs, and their 
evolution both with complex parameters of the Hamil- 
tonian and with a non-integrable complex perturbation. 
On this basis, we prove that integrable models with at 
least two parameter-dependent integrals of motion (IMs) 
have a level crossing in the complex-extended parameter 
space, providing a clear signal of their integrability. Fur- 
thermore, the inclusion of a non-integrable perturbation 
splits each level crossing into two EPs transforming dra- 
matically the topology of the Hilbert space close to the 
level crossing. 

Let us begin by briefly reviewing the RG models 
which are based on the SU(2) algebra with elements 
Ki, , and Kf, fulfilling the commutation relations: 
[K+,K„] = 6u>K° , [K", Kp] = ±5 W K± . The indices 
I, I' refer to a particular copy from a set of L, SU (2) alge- 
bras. Each SU (2) algebra possesses one quantum degree 
of freedom. Therefore, a quantum integrable model re- 
quires the existence of L independent, global operators 
that commute with each other. These operators, which 
need not be hermitian, are the IMs. In the following, we 
will work with the rational family of RG models whose 
IMs are 
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where g and £; are L + 1 arbiratry parameters. The IMs 
(P) satisfy [Ri, Rj] = for all pairs i, j. 

There is a profound relation between quantum integra- 
bility in finite systems and the existence of level cross- 
ings. In the context of the 6-sites Hubbard model, this 



problem has been addressed by Yuzbashyan et al. 15 1 
who showed that a level crossing implies the existence 
of two independent parameter-dependent IMs. Here we 
complete this analysis by showing that a system with 
at least two parameter-dependent IMs has of necessity a 
level crossings in the complex plane. 

Let us assume a Hamiltonian H (g) depending linearly 
on a parameter g. H(g) itself is the parameter-dependent 
IM. If g is complex then the Hamiltonian and eventually 
other parameter-dependent IMs will be non-hermitian. 
The parameter-dependent IM Q(g) commuting with the 
Hamiltonian: [H(g),Q(g)] = 0, will be independent of 
the Hamiltonian if it cannot be expressed as an entire 
function of H(g) and g, i.e., Q(g) ^ f(g,H(g)). Let 
us assume that n is the dimension of the Hilbert space 
and Ex(g), ■ ■ ■ , E n (g) and qi(g), ■■ ■ , q„(g) are the corre- 
sponding eigenvalues in the basis in which both operators 
are diagonal. If Q(g) is an entire function of H(g), then 
it can be expanded for any complex g value as: 
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This set of equations has always a solution unless for 
some value g — go a pair of equations {k, k'} have 
Ek(go) = Ek>(go), but q k (ga) ^ qk'(go), implying a dou- 
ble degeneracy in the Hamiltonian but not in the second 
IM. 

The minimal rational RG model {1} allowing for level 
crossings should have at least three SU(2) copies. The 
reason is that the sum of the IMs |T]) is a global conserved 
symmetry: A" = X^iLi K?> commuting with all IMs and 
independent of the parameter g. Hence, we are left with 
two parameter-dependent IMs and, as shown above, at 
least two parameter-dependent IMs are required for hav- 
ing level crossings. 

In what follows, we will use the pair representation of 
the SU(2) algebra leading to pairing Hamiltonians. The 
elementary operators in this representation are the num- 
ber operators Nj = J2m a ]m a jm an d the pair operators 

Aj = Em a jm a jm' where j is the total angular momen- 
tum and m is the z-projection. The state jm is the time 
reversal of jm. 

The relation between the operators of the pair algebra 



and the generators of the SU(2) algebra is: K$ 

jttj , Kf = (Kj^y = where Qj is the particle 

degeneracy of level j. With this correspondence, one can 



introduce an integrable 3-level pairing Hamiltonian as: 
H(g) = 2j2eiRi(9) + C = ]>>^ + <?]T A* A,, (3) 



where C is an irrelevant constant and £j (i = 1,2,3) 
are the single-particle energies. One can see that H(g) 
itself is a parameter-dependent IM. As discussed above, 
the sum of the IMs (Q} yields a parameter-independent 
IM, the particle number: N = 2 £ ; R t (g) + i J2i n Z ■ The 
second parameter-dependent IM can be chosen as linearly 
independent from the two other IMs. The simplest choice 
is Ai. If £\ = 0, the second parameter-dependent IM 
becomes: 



Q(g) = 



l 



'fi 2 , Q 3 
5 1 — + — 

£2 £3 



9 \i 

1 

£.3 



\(aW 

\ (a{A 3 + AlAi 



Ni 
2 

AAi 



A«2 gfti N 3 
£ 2 2 e 3 2 

+ NxN 2 



NiN 3 



(4) 



The position of all degeneracies in the complex-g plane 
are indicated by the roots of the coupled equations: 

dct [H (g) - EI] = ; A de t [H (g) - EI] = 0. (5) 

By eliminating E from these two equations, we are left 
with the discriminant D(g), a polynomial in g of degree 
M < n[n — 1). The discriminant can be written as [9(: 



D(g)= J] [E m (g)-E m >(gf 
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where E m (g), E m r(g) denote the complex eigenvalues of 
H(g). The eigenvalue degeneracies E m (g) — E m /(g) at 
g = g a {a = 1,...,M), can be found numerically by 
looking for sharp minima of D(g). It turns out that for 
real values of Sj, the roots are real or complex conjugate 
pairs. One finds two kinds of solutions in quantum inte- 
grable models: (i) single roots corresponding to EPs that 
are common to all IMs, and (ii) double- (multiple-) root 
solutions which indicate non-singular sharp crossings of 
two (or more) levels with two (or more) orthogonal wave 
functions, related to the existence of at least two inde- 
pendent IMs. 

Fig. [T] shows an evolution of two level crossings, the 
double-roots of D(g), in the complex-g plane as a func- 
tion of the energy of the third single-particle level £3 
(£3 > £2). The systems has 4 pairs of fermions in a va- 
lence space with level degeneracies f2i = 6, 02 = 4, = 
2 and e\ = , £2 = 1. In the limit £3 — * 00, this sys- 
tem decouples effectively into the two 2-level models: the 
first one with level £3 occupied and the second one with 
£3 empty. In this limit level crossings are forbidden and, 
indeed, the two level crossings that appear for finite val- 
ues of £3 move to ±00. With decreasing £3, two level 
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FIG. 1: Collision and subsequent scattering of two level cross- 
ings as a function of the energy £3 in a complex-extended in- 
tegrable 3- level pairing Hamiltonian ©. Points are plotted in 
a descending order of £3. For more details, see the discussion 
in the text. 



crossings coming from ±00 approach each other in the 
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real axis up to a critical value £3 = 1.8499 where they 
coalesce. The level crossing at this point corresponds to 
the quadruple-root of D{g). For £3 < £3 this cross- 
ing splits again into the two double-root level crossings 
which move into the complex- plane. The presence of 
such level crossings in the complex plane is a clear signa- 
ture of quantum integrability, and shows the necessity of 
extending the demonstration of operator independency 
given in @ to the whole complex parameter space. 

EPs associated to single roots of the discriminant D(g), 
unlike level crossings, are common to all parameter- 
dependent IMs including the Hamiltonian. It should be 
also noted that EPs appear in the quantum integrable 
model even though no manifestation of level repulsion is 
expected and the spectral fluctuations of the hermitian 
Hamiltonian obey a Poisson distribution [lij]. In that 
sense, level repulsion may be a sufficient but not a nec- 
essary condition for the appearance of EPs. 

Another feature that we found in this 3-level pairing 
model, as well as in more general multi-level pairing mod- 
els, is the reduction of the total number of discriminant 
roots whose maximum value is M max — n{n — 1). In the 
particular case shown in Fig. 1, M max — 20 but we found 
16 roots consisting of 2 level crossings (double roots) and 
12 EPs (single roots). 

In the following, we generalize the 3-level pairing 
Hamiltonian ^ to study effects of non-integrability: 

H{g) = E + (9 E ~ (1 - 09 E N i ■ ( 7 ) 

i ij i 

For £ = 1, Eq. Q corresponds to the integrable Hamil- 
tonian ([3]). For £ = 0, the Hamiltonian (J7JI is also inte- 
grable with the number operators playing the role of 



parameter-independent IMs. In the interval: < C < 1, 
the Hamiltonian ([7]) is non-integrable, i.e. it does not 
possess an independent IM other than the Hamitonian 
and the total number operator. Two main features char- 
acterize the emergence of non-integrability, firstly the 
level crossings break into pairs of EPs and secondly the 
missing roots of the discriminant come into play as EPs 
from 00. 




FIG. 2: The evolution of level crossings and EPs of the 
complex-extended 3-level pairing Hamiltonian © in the in- 
terval < £ < 1. Circles and squares denote the position of 
degeneracies at £ = 1 and 0, respectively. The double circles 
depict the double-root level crossing at £ = 1. The triple cir- 
cle shows the quadruple-root level crossing corresponding to 
the coalescence of two double-root level crossings. For more 
details, see the caption of Fig. [T] and the discussion in the 
text. 



Fig. [2] shows the global pattern of level crossings and 
EPs as a function of the parameter ( (0 < ( < 1) for 



(cr) 



the 3-level system of Fig. 1 in three regimes: £3 > e 

(£3 = 7/3), £3 = £ 3 cr) , and £ 3 < £ 3 cr) (£3 = 3/2). Ener- 
gies of levels T' and '2' are fixed at: e\ = 0, e% = 1. Only 
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the lower half-plane of g is shown where all eigenvalues 
are either discrete states on the real-g axis or else decay- 
ing resonances. Complex conjugate degeneracies situated 
in the upper half-plane (Tm(g) > 0) correspond to cap- 
turing resonances. In the limit of £ = 1, the two level 
crossings for e 3 = 7/3 are located along the real-g axis. 
One may also notice a quadruple-root level crossing on 
the real-g axis at £ 3 = £g Cr ' (see also Fig. [T]). Moreover, 
one can see the location of 6 EPs of the integrable pairing 
Hamiltonian ^ in all regimes of £3. 

With decreasing £, one observes several distinct effects. 
Firstly, each double-root level crossing at £ = 1 breaks 
into a pair of EPs. For £3 = 3/2, two EPs resulting from 
this fragmentation follow independent trajectories in the 
complex-^ plane and end up in two different level cross- 
ings at C = 0. For £3 = 7/3, the level crossing at £ = 1 
breaks into two EPs symmetrically with respect to the 
real-g axis. Since this symmetry is conserved for all £, 
these EPs end up in the same level crossing for £ = 0. 
Secondly, roots that are missing in the integrable limit 
(C = 1) appear from g — 00 and end up in level crossings 
for C -> 0. Thirdly, in the limit £ ->• 0, all EPs either 
collapse in different level crossings at the real-g axis or 
escape to infinity. The double level crossings found in 
this limit correspond to the two different eigenvalue de- 
generacies. For £ 3 = 3/2, one can see three EPs converg- 
ing on each side of the real-g axis to a single point. At 
this sextuple-root of the discriminant, one finds a sharp 
crossing of three eigenvalues. 

The degree of the discriminant, which for ( = 0,1 
equals 16 in all regimes of £3, becomes M = n(n— 1) = 20 
in the non-integrable regime. Notice the absence of EPs 
in the integrable case C = reflecting the fact that the 
Hamiltonian is diagonal in the original basis for any g 
value. 

In conclusion, we have shown that a system with two 
independent parameter-dependent IMs has at least one 
level crossing in the complex parameter space. We have 
used a minimal integrable pairing model consisting of 
three levels to exemplify this issue. Moreover, we have 
found that the degree of the discriminant is reduced in 
this integrable limits, giving rise to a lower number of 
EPs, but still a fraction of them persists in spite of the 
fact that there is no manifestation of level repulsion. If 
integrability is broken by the inclusion of a non-integrable 
perturbation, all level crossings split into EPs and other 
EPs come into the complex plane from 00 recovering the 
maximal degree M max = n(n— 1) of the discriminant. We 



conjecture that these two unexpected properties, level 
crossings in the complex-g plane and the reduction in 
the number of EPs, uniquely define a quatum integrable 
system. If this conjecture proves to be true, then it might 
be particularly useful in studies of finite systems with 
small dimensional Hilbert spaces, where the analysis of 
spectral fluctuations is unreliable. More work has to be 
done to elucidate the relation between generic integrable 
Hamiltonians and the missing roots of the discriminant 
in order to confirm this conjecture. 
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